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Abstract 

We propose an interpretation of the random close packing 
of granular materials as a phase transition, and discuss the 
possibility of experimental verification. 
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Introduction 

The phenomenon of random close packing was popularized by Bernal [B] as 
part of his effort to model liquids. The experiments which clarified many of the 
issues were performed by Scott et al [S,SC,SK], and showed the following. If a 
large number of monodisperse hard spheres, for instance steel ball bearings, are 
gently poured into a container, the volume fraction will be roughly 0.61. If the 
container is repeatedly shaken vertically, this density rises to about 0.64, and care- 
ful protocols lead to reproducible lower and upper limits on the volume fraction, 
called, respectively, random loose packing (0.608 ± .006) and random close packing 
(0.6366 ±0.0005) [SK]. Scott et al [SC] noted that volume fractions beyond 0.64 (up 
to approximately 0.66) could be obtained if the material was cyclically sheared, and 
that this result was accompanied by small crystal-like clusters of spheres. This was 
confirmed and explored by Pouliquen et al [ND] (who also performed experiments 
employing horizontal shaking [PN]), obtaining volume fractions up to 0.70, again 
accompanied by crystal-like clusters of spheres. (Recall that the densest possible 
packing of monodisperse spheres has volume fraction tt/vTS ~ 0.74.) 

In light of the above, the volume fraction 0.64 is generally described as the 
boundary between two regimes for granular matter: at volume fractions below 0.64 
the structure of the material is random, while above 0.64 it has some order, as 
represented by the crystal-like clusters which appear. 

This general understanding has recently been questioned in an influential paper 
by Torquato et al [TT] , in which the authors claim to "have shown that the notion 
of RCP (random close packing) is not well defined mathematically". In contrast 
we propose an unambiguous meaning to such a boundary between disordered and 
ordered states of granular matter, as a boundary between well defined phases, to- 
gether with a mathematical model of traditional form and an experimental test of 
our interpretation. 

Analysis 

Our proposed characterization of random close packing is motivated by proper- 
ties of the hard sphere model of classical statistical mechanics. In that model point 
particles, with the usual position and momentum degrees of freedom, interact only 
through a hard core: no pair may approach closer than some fixed separation D, 
and the particles evolve dynamically through elastic collisions of imaginary spheres 
of diameter D surrounding their positions. Our interest in the hard sphere model 
stems from the demonstration by Alder and Wainright [AW] , by molecular dynamics 
simulations, that the model exhibits a first order phase transition between a fluid, 
which exists at volume fractions below 0.494 ± 0.002, and a solid, believed to be 
crystalline, which exists at volume fractions above 0.545 ±0.002 [HR]. Between 0.49 
and 0.54 there is a mixed phase. Using the canonical ensemble we integrate out the 
momentum variables and consider the "reduced" probability distributions on the 
phase space of the position variables alone, in the inflnite volume limit. They are 
the inflnite volume limits of the uniform distributions on packings for flxed volume 
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fraction. In particular the distribution p(m) / of the mixed phase at volume fraction 
/, 0.49 < / < 0.54, is represented by an average of the distributions of the pure 
phases: p{m) f = cpo.54 + (1 ~ c)Poa9i where po.49 is the distribution of the highest 
density fluid, po.54 is the distribution of the lowest density solid, and < c < 1 
is such as to produce the volume fraction /, namely c = (/ — 0.49)/0.05. (This 
is merely a statement of the fact that distinct phases separate when coexisting in 
equilibrium, each occupying a well defined volume [LL].) Note that 0.49 is the vol- 
ume fraction of "freezing" . Therefore assuming, as generally believed, that the solid 
phase is crystalline (in fact face centered cubic [W,BF]), we interpret 0.49 as the 
"highest random density" among monodisperse spheres. (These structural features 
have been confirmed not only by many computer simulations but also in experi- 
ments with appropriate colloids [RD].) Intuitively, at any volume fraction above 
the freezing point there is a nonzero probability of (seeing) an infinite, ordered 
crystal. It is the use of the infinite volume limit, together with the probabilistic 
formalism, which produces a sharp phase transition between disorder and order in 
equilibrium statistical mechanics [FR]. 

We emphasize that there are sphere packings with packing fraction d > 0.49 
which might well be described as random, for instance packings corresponding to 
any metastable extension of the fluid branch of an isotherm. However the total 
of all such packings has probability zero with respect to the (infinite volume limit 
of the) uniform distribution on packings with packing fraction d. Our use of this 
distribution as the touchstone of relevancy is in accord with its common appearance 
in statistical physics and probability theory; the best justification one can give is 
that it is commonly found that practical sampling of phase space seems to occur in 
this way, in particular in the natural dynamics of matter in thermal equilibrium. 

In summary, the hard sphere model, and its physical realization in colloids, 
exhibits the basic ingredients needed to make sense of the granular phenomenon of 
random close packing: the volume fraction 0.49 separates the fiuid phase of random 
packings from the mixed phase in which crystalline order begins to appear. While 
this is not directly applicable to the granular matter which is our proper subject, it 
nonetheless shows that the intuitive notion of random close packing is not inherently 
inconsistent as claimed in [TT]. 

We now turn to granular matter. The traditional hard sphere model does 
not include the effects of gravity and cannot represent the properties of granular 
matter, in particular that of random close packing. However a slightly modified 
ensemble framework has been proposed as a model for granular matter. Specifically, 
in the original proposal of Edwards ct al [EO] one uses a uniform distribution 
on those static monodisperse sphere packings, of fixed volume fraction, which are 
mechanically stable under gravity. One can add friction to the spheres and perhaps 
other restrictions besides volume fraction; adding a condition of fixed pressure might 
be useful, though it is not clear if pressure is isotropic in granular materials. 

As is true for solids in thermal equilibrium, in order that an ensemble method 
be appropriate for {nonequilibrium) granular matter it is important to restrict the 
protocols used to produce beds of granular matter at fixed volume fraction (and 
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perhaps pressure etc.) One feature that is necessary is that the protocol be "his- 
tory independent" [NK] in that it give equivalent results starting with beds orig- 
inally prepared in any manner. There are three types of protocols which claim 
to produce history independent beds of monodisperse granules with well defined 
volume fractions: vertical vibration (often called "tapping"); fluidization followed 
by sedimentation; and cyclic shearing. These methods have produced history in- 
dependent beds with volume fractions in the following ranges: 0.605 to 0.625 by 
vertical vibration [NK,RR]; 0.685 to 0.70 by cyclic shearing [ND]; and 0.57 to 0.62 
by fiuidization/sedimentation [SG,SN]. 

It is known from Schroter et al [SN] that granular beds prepared by fiuidiza- 
tion/sedimentation undergo a phase transition, as volume fraction is varied, at ap- 
proximately 0.60 volume fraction, as measured by two different responses to shear. 
Given the mathematical similarity between the hard sphere model and the granular 
model of Edwards on the one hand, and the experimental similarity between the 
sharp freezing transition in colloids [RD] and the abrupt appearance of crystalline 
clusters in the experiments of Scott et al [SC] and of Pouliquen et al [ND] on the 
other hand, we predict that history independent granular beds would show another 
phase transition: a first order phase transition, with a mixed phase for volume 
fractions between 0.64 and 0.74, again exhibited through the response to shear or 
other mechanical probe. Analogously to the hard sphere model, the distribution 
for granular matter in the mixed phase would be a mixture, with one component, 
at volume fraction 0.74, representing a crystal, and the other, at volume fraction 
0.64, representing a disordered phase. This would give a well defined meaning to the 
phenomenon of random close packing of granular matter just as the freezing density 
defines a similar concept for the hard sphere model or for hard sphere colloids. 

We note the connection between this proposal and that of Kamien and Liu 
[KL], which also uses the hard sphere model to understand random close packing. 
In [KL] random close packing is associated with the end point of a metastable 
branch in the hard sphere phase diagram, while we use the hard sphere model 
only to predict behavior in a related but different ensemble, of packings which are 
mechanically stable under gravity, and in particular we predict a phase transition 
at volume fraction 0.64. 

History independent experimental protocols have not yet produced beds with 
volume fraction in any interval containing the volume fraction of interest, 0.64; 
this will be necessary before our prediction can be checked against the behavior 
of granular matter. Alternatively it might be possible to test the prediction by 
realistic, history independent computer simulations. However it should be noted 
that fifty years of computer simulations of the hard sphere model have not yet 
been able to demonstrate the appearance of crystals at its freezing transition, so 
one should not be too optimistic that an order/disorder transition could be seen in 
simulations of granular models. 
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